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A variational expression for the admittance of a hollow cylindrical antenna IV<[ at the 
center has been obtained. Contrary to normal practice, the integral equation and the re- 
sulting variational expression have as the unknown function, the tangential electric field in 
the tubular region extending from the antenna ends to ± °° . A plausible trial field is chosen 
and the relevant integrals are evaluated asymptotically for long, thin antennas. Some 
numerical results are given which show reasonable agreement with the experimental values 
of [izuka et al. [1963], and also with the theoretical results of W'u [1901]. 



1. Introduction 

Recently there has been considerable interest 
shown in the admit t mice of long cylindrical antennas. 
Vainshtein [10591, by assuming that the current 
consisted of incoming and outgoing waves with ampli- 
tudes which were slowly varying functions of posi- 
tion, reduced the problem to the solution of Volterra 
integral equations which could be solved exactly. 
Another approach was used by Wu [1961] who 
approximated a long antenna by a semiinfinite one. 
The latter problem was then solved by Wiener-Hopf 
techniques. Extensive tables of numerical values 
were presented by Wu. Quite a different method 
was used by Chen and Keller [1962], who assumed 
the current to be composed of a wave emanating 
from the feed point plus waves reflected from the 
ends. The current reflection coefficient at the end 
of a tube was known from previous work. With a 
knowledge of this quantity the authors were able to 
obtain the admittance without formally solving an 
integral equation. Essentially the same method 
was used earlier by Hallen [1956, 1962]. Studies on 
the infinite antenna were carried out by Hallen 
[1948], Papas [1948], Duncan [1962], and Kunz [1963]. 
Measurements on long antennas have been made by 
Iizuka et al. [1963], whose results agree quite well 
with the theoretical values of Wu [1961]. Vari- 
ational net hods were applied to the antenna prob- 
lem by Storer [1950] and by Tai [1950], who obtained 
variational expressions for the impedance in terms 
of the unknown current distribution on the antenna. 
Such a formulation is expected to be most accurate 
for short- antennas. To obtain an integral equation, 
and from it a variational expression whose approxi- 
mate solution is most accurate for long antennas, 
one should use the tangential electric field in the 



"aperture" as the unknown function. 1 Here the 
aperture is understood to be the tubular extension 
reaching to ± oo from each end of the antenna. To 
the writer's knowledge an integral equation contain- 
ing the electric field as unknown lias never been 
used, and this approach forms the substance of this 
paper. 



2. Integral Equation and Variational 
Expression 

The antenna is assumed to be a thin-walled, 
perfectly conducting tube occupying the region 
P=a, —h<z<h, of the (p, <t>, z) cylindrical coor- 
dinate system. At 2=0 the antenna is driven by a 
voltage V applied across an infinitesimal circum- 
ferential gap. The only electric field component in 
the gap is given by 



E z =-Vb{z)e i<at . 



(1) 



The space surrounding the antenna is divided into 
two regions, denoted "1" for p<^a and "2" for p>tf. 
Field components in these regions will henceforth 
be distinguished by appropriate superscripts. 

It is assumed that all field components are deriva- 
ble from the ^-components of electric hertz vectors. 
In either region, after dropping the time-dependence 
factor, 



r* co 



A^ZtifrpY^da, 



= 1,2; 



(2) 



1 This seems a reasonable extension of a statement of Levine [1955] p. 38, who 
was concerned with diffraction by a circular aperture. 1 am indebted to R. F. 
Millar for supplying this reference, and for a helpful discussion of this point. 
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where /3 2 =k 2 —a 2 , the A^ot) are amplitude functions, 
and the Z t are Bessel functions appropriate to the 
region. The contour of integration is assumed to be 
indented downward around the branch point a =—k, 
and upward at a=k. From (2), the tangential elec- 
tric and magnetic field components are easily 
obtained: 



J 00 
— 00 

?=— twof pA i {a)Z' i (fip)e iaZ da, 

*7 — oo 



(3) 



(4) 



the dash on the Bessel function indicating derivative 
with respect to the argument. At p—a, E z must be 
continuous for all z, thus 



A 1 (a)Z 1 (fia)=A 2 (a)Z 2 (Pa). 



(5) 



Also E 2 =0 for p=a, — h<z<h, with the exception 
of 2=0; so that, upon taking the inverse Fourier 
transform of (3) and using (1), 



2T(3 2 A i (a)Z i ((3a)= f e 

J \z\>h 



e- iaZ E z (z)dz-V. (6) 



In the region z^>h, p=a, H^ must be continuous, 
hence from (4) 



r 



Pe^ [A l Z[{fia)-A 2 Z f 2 {fia)] da=0, \z\>h. (7) 



For \z\<Ji, p=a the total current I t on the antenna 
is given by I t =—2Ta[H£ ) —H£ ) ], which from (4) is 

I t (z)=2irio)e a f ^^^(^Z^^-A^Z^a^da, 

\B\<h. (8) 

Use is made now of (5) and (6) to remove Ai(a) and 
A 2 (a) from (7) and (8) with the result 



j. 



e iaZ da 



_„/3 2 Jo03a)tff(/3a) 



if ■ 



e-^'E&'W-V y=0, \z\>h (9) 



? 



, -"-/«(s) / 

2coe 

M<h; 



(10) 



where Z x has been replaced by J and Z 2 by H^ 2) , 
and the Wronskian relation J H^ 2)/ — H^ 2) Jq= 
—2i/(irl3a) has been used. Interchanging the order 
of integration and writing E(z) for E z (z) yields 
the integral equation for E(z) 

f E(z')K(z-z')dz'=VK(z), \z\>h, (11) 

J\z'\>h 



and the relation 



f E{z')K{z-z')dz> f =VK(z)-^ 

J\z'\>h ^0. 



(2) 

7 

coe 



where 



/» oo 



e iaZ dc 



P 2 J (t3a)H™(Pa) 



\<h; 
(12) 

(13) 



To obtain a variational expression, (11) is first 
multiplied by E(z) and integrated with respect to z 
over \z\^>h, to give 

(( E(z)E(z')K(z-z')dzdz' 



=V f E(z)K(z)dz. (14) 

J\z\>h 



If z is allowed to approach zero in (12) one obtains, 
since K( — z)=K(z), 



i 



E(z)K(z)dz=VK(0)- 7 ^^' (15) 

\z\>h ^we 



Equation (14) now is divided by [fE(z)K(z)dz\ 2 
and (15) used to introduce I t (0), giving 



2cO€ T' 7 



■K(oy- 



[ f E(z)K(z)dzl 

LJ\z\>h J 



(( E{z)E{z')K(z-z')dzdz' 

J J\zl\z'\>h 

(16) 



This is a variational expression for the quantity 
on the left. The proof of this is omitted, but 
follows the standard procedure. An interpretation 
of the quantity K(0) is now required. This is 
readily obtained from (12). If A-^>°o the integral on 
the left tends to zero, so that 



K(z) = icI?{z)l2ue*V, 



(17) 



where /," is the total current on an infinite antenna. 
If z approaches zero, then 



K(0) = irY?/2ue , 



(18) 



with Y? the total admittance at the feed point of an 
infinite antenna. Then (16) becomes 



Y t -Y? 



[fE(z)K(z)dz] 



2o)6 

"V 7 fE(z)E(z')K{z-z')dzdz'~ 



(19) 



This is the variational expression in its final form. 
It relates the admittance of a finite but long antenna 
to that of the infinite one. The latter has been the 
subject of considerable research, to which Papas 
[1948], Hallen [1948], Duncan [1962], and Kunz [1963] 
have all contributed. 
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A remark about (19) is in order. Although the ad- 
mittances given are total ones; that is, they include the 
effects of currents inside the tube, the expression (19) 
is approximately true also for the external admit- 
tances (called Y and Y m ), when the tube diameter is 
reasonably small. This follows because the interior 
currents are excited by evanescent waveguide modes 
which are rapidly attenuated by distance, and which 
consequently are hardly affected by the antenna 
ends. Thus the interior currents are almost equal for 
infinite and long antennas and hence cancel in (19). 

3. Evaluation of the Variational Expression 

The problem now arises of choosing a suitable field 
distribution E{z) to use in (19). There are two ob- 
vious requirements which this function should satisfy: 



it must behave like 



-ik\z\ 



l\ 



for laro-e \z\. and it 



should vary as [\z\ — h] l/f * near \z\ = h. Both can be 
achieved by taking 



E(2)=\z\-\z 2 -h 2 )- 1/2 cxp (-ik\z\). 



(20) 



A constant multiplier might have been included in 
(20), but was omitted because the variational expres- 
sion is independent of this term. 

When (20) is inserted in (19), integrals of the fol- 
io win % kinds arise: 



7 _rr cx\)[-ik(\z '\)]K{z- z r )ihdz\ 

JJ\z\,\z'\>h 



I 



=L 



(2 2 -/r) 1/2 (s -h 2 ) ll2 \zz'\ 
exp (— ik\z\)K{z)dz 

>h {z 2 -h 2 )" 2 \z\ 



(21) 



(22) 



The approximate evaluation of these integrals is car- 
ried out in the appendix, under the assumptions that 
the radius is small in terms of wavelength and that 
the length of the antenna is great. The results are: 



i^ 



where 



lit" 



vir 



2k' 2 fi r L 2 exp t -4 *^ +FFZ exp (~ 2ikh "> > ( 23 ) 
/^(^/kWy'Lr 1 exp (-2ikh-iT/4), (24) 

Z,„=ln (iirka) 
Z,=ln (-2ih/rka 2 ) 
L 2 =ln(-4ih/Tka 2 ), 



and r= 1.781. . . . If the results (23) and (24) are used 
in (19), the following expression for the admittance 
Y—Y" is obtained 



Y-Y"= 



-4irY n L 2 



mi-(2L s /vL ) exp (2ikh)] 



(25) 



where F =(120x) ! mhos is the admittance of free 
space. 



Expression (25) closely resembles, in form, HilleV 
result [1962, p. 485]. Apart from a difference in the 
argument in one of the logarithms, which is of the 
same magnitude as higher order terms omitted in the 
expansion of Y—Y m , the only difference is the re- 
placement of a factor 2/tt in the denominator of (25) 
by a factor § in Hallen's work. This difference has 
quite a large effect near resonance. 



4. Discussion of Results 

To compare (25) with the results of other workers, 
an expression for Y 00 is needed. A difficulty is that 
F 00 turns out to be infinite, as noted by Infeld [1947] 
and others. This point has been discussed by Dun- 
can [1962], who proposes a "smoothed" form of the 
admittance for an infinite antenna. In the present 
paper values interpolated from Duncan's table 3 are 
used. Two different radii are employed : a/\=0.0085 
and a/\=0.001191. The values of Y m from Dun- 
can's paper corresponding to these were found to be 
0.00300 +0.001 87i and o'.00200+0.00084i mhos, re- 
spectively. 

The results of the calculation based on (25) and 
carried out for three different ranges of antenna 
lengths are shown in figures 1, 2, and 3. Some ex- 
perimental data, for long antennas have been given 
by fizuka et al. [1963], but owing to the sniallness of 
their graphs, it is not easy to compare their results 
directly with calculations based on (25). However, 
their curves follow closely Wu's [1961] theoretical 
values, which arc tabulated, and which have been 
included in figures 1, 2, and 3. It will be observed 
that the present results for the range 51.0<M<52.6 
(fig. 3) are displaced to the right of Wu's, which in 
turn are on the right of the experimental curves of 
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Ficjure 1. Admittance of a cylindrical antenna, &/\ = 0.00119l. 
1.0 <kh< 2.0. 
Dots indicate Wu's theoretical values. 
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5. Conclusion 

In conclusion, a simple variational solution for the 
admittance of a cylindrical antenna lias been ob- 
tained, involving the longitudinal electric field, in the 
extension of the antenna as unknown. This formula- 
tion is most accurate for long antennas, in contrast 
with the variational method in terms of current, 
which is more suitable for short antennas. For long 
antennas it is easier to choose a plausible electric 
field than a current, since the latter involves waves in 
two directions with an unknown reflection coefficient . 
Using a reasonable trial field, the first term of an 
asymptotic series for the admittance is obtained. 
This gives fairly good agreement with experiment 
and with Wu's theory [1961]. More accuracy will 
be obtained if extra terms in the asymptotic develop- 
ment are included, provided the antenna is suffi- 
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The author thanks Mrs. M. Steen for carefully 
doing the calculations, and E. V. Jull for reading and 
criticising the manuscript. 

6. Appendix. Evaluation of the Integrals 

The integral I x of (21) is first considered. By 
halving the ranges of integration and taking account 
of the evenness of K(z) , I x can be written in the form 

1 -r da 
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1 "Jo 0>Jo(pam(f3a) 

XL z(z 2 -h 2 ) y ' J 

The term in braces can be expressed in terms of 
Hankel functions through the relation 

giving approximately, for large h 

+HV[(k-a)h}} 2 - (27) 

For reasonably large h, H [(a-\-k)h] can be replaced 
by its asymptotic value. It is desirable to approxi- 
mate the other Hankel function the same way, but 
it is not immediately clear that this can be done, 
since its argument is small near a=k. However, 
the contour can be deformed into the first quadrant 
away from a — k without crossing any singularities, 2 
hence on the new contour the approximation is 
valid if h is sufficiently large. The final step con- 
sists of deforming back to the original contour with 

2 The poles of the integrand are the zeros of Jo((3a) which all lie on the contour 
Im (0) =0. By virtue of the time dependence chosen these lie in the 2d and 4th 
quadrants. 
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the result 

1 h? e Jo pJMa)mm\k+aT r k-a r flJ 

- (-irtA ,, )«- ,a *l/u+/u+/is]. (28) 

Integrals of the type I u take their values mainly 
from the neighbourhood of <* = &; so that 



T »* P 



1 e- 2tah da 
/3 2 ln(ir^ffl)' 



(29) 



where the Bessel functions have been replaced by 
their small argument approximations. The integral 
in (29) has for its leading term 



In**. 



2F Lo 



(30) 



[Hallen 1948]. 

The integral 7 12 , though of similar appearance, 
takes quite a different form. Since e 2ia *-»0 for 
Im (a)^oo, the contour can be moved one along the 
positive imaginary axis: 



/i 



,-i J. (*= 



e- 2 ^f 



ir)(^+f 2 )Jo[a(^+f 2 ) H ]H^[a(A 2 +r 2 ) M ] 

(31) 



Watson's lemma (see Jeffreys and Jeffreys [1962], 
p. 501) can now be applied, with the result I 12 = 
l)(l/h), so that 7 12 , to the first approximation, is 
negligible compared with I v . Strictly speaking, the 
poles which lie near, and in the absence of loss on, 
the path of integration should be accounted for, 
but if a is small, they occur for large values of f and 
thus contribute very little to 7 12 . 

There remains I n . The contour of integration is 
first changed to the positive imaginary axis, with the 
result 



--2i 



/»oc 

Jo 



di 



(F+r 2 ) ?4 e7 [a(* 2 +f 2 ) K ]^'[a(A 2 +f 2 )^] 



(32) 



The range of integration is now split at f =&, so that 
in (0,6) the small argument approximation to the 
Bessel function applies. The remaining integral on 
(b, oo ) y being of order 1/b 2 and hence 0(# 2 ), is ne- 
glected. Then 



'**-*£ 



<n 



(A 2 +f 2 ) 3 =ln[|ira(F+f-')'fl 



(33) 



In (33) the logarithmic term is a slowly varying 
function, and without large error f in it can be re- 
placed by a constant lying in the range (0, b) . The 
value chosen is not very important since any error 
is of the order of the terms omitted in the expansion. 



It is reasonable to take f=0, however, since the term 
(& 2 +f 2 ) is smallest there. The integration remaining 
is elementary, giving 

/ 13 ~-7r6/t 2 i (* 2 +6 2 ) 1/2 . 
Recalling that b is large, approximately 

Iu=-t/Jc 2 L . (34) 

When the terms are gathered together, the final 
result is 



/:~- 



kWL n 



(35) 



The integral L given in (22) is now evaluated. 
Following the course outlined for I u the integral I 2 
becomes 

+H?[(k-a)h]} (3G) 

which, as before, can be further simplified to 

(la 



I 



m^L 



C e ia " e- iah ~\ 



(37) 



The terms in (37) are similar to /,, and I l2 in (28), and 
may be handled in like manner, with the result 

I 2 <^(t 5 /JcW) h U ] exp (-2ikh-i<ir/4) . (38) 
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